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Abstract 

Abstract: In the area of bipedal locomotion, the spring loaded inverted pendulum (SLIP) model 
has been proposed as a unified framework to explain the dynamics of a wide variety of gaits. In this 
paper, we present a novel analysis of the mathematical model and its dynamical properties. We 
use the perspective of hybrid dynamical systems to study the dynamics and define concepts such 
as partial stability and viability. With this approach, on the one hand, we identified stable and 
unstable regions of locomotion. On the other hand, we found ways to exploit the unstable regions 
of locomotion to induce gait transitions at a constant energy regime. Additionally, we show that 
simple non-constant angle of attack control policies can render the system almost always stable. 
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I. INTRODUCTION 



One of the most accepted mathematical models for bipedal running is the spring loaded 
inverted pendulum (SLIP, for an extensive review see[l(). In a similar fashion, the rigid 
inverted pendulum has been extensively used to model bipedal walking^. In 2006, Geyer 



et al.[3| propose the SLIP model as a unifying framework to describe walking as well as 
running. The unified perspective proves useful for accurately explaining data from human 
locomotion jij]. Additionally, it allows describing both gaits (walking and running) in terms 
of dynamical entities observed in a discrete map, obtained by intersecting the trajectories 
of the system with a predefined section of lower dimension. Geyer associates these entities 
with limit cycles of the hybrid dynamical system^, Q and named their attracting behavior 
as self-stabilization. Though the nature of the observed dynamical properties is not yet 
clarified, those results emphasize that bipedal locomotion may be dictated solely by the 
mechanics of the system. As a consequence, the control necessary for locomotion is thus 
reduced to the swing phase of the leg, showed in Fig. [1] between points A and B. The most 
popular control policy is to produce touchdowns at constant angle of attack a (CAAP (a)), 
i.e. the angle spanned by the landing leg and the horizontal. 

In the last decade, many energy-efficient bipedal walking machines have been developed. 
Through careful design, they exploit the passive dynamics of their own body to move forward, 
requiring little control or none|6NlO|. However, the construction of bipedal machines capable 
of exploiting passive dynamics in different gaits remains an unsolved engineering challenge. 
In this context, Geyer et al. [3[ report that, in the SLIP model, it is not possible to have 
multiple gaits at the same energy. The results are based on simulations that do not cover all 



possible initial conditions of the system. In addition, Rummel et al. Uj prove that walking 
and running is possible at the same energy level. They use a new map that allows comparing 
different gaits with ease. The map is defined at the vertical plane crossing the landing point 
of the foot (Fig. [1]). In this way, they find the self-stable regions, but their intersection is 
empty. To concretize these ideas, let us describe this region for the running map 1Z. 

E* = {x\ x eSA(3a\ x = ll a (x))}, (1) 

where the subscript in TZ a denotes running using CAAP (a) and S denotes the section where 
the map is defined. Therefore, if for different gaits these stable regions do not intersect, e.g. 



fl E^ = 0, we conclude that a transition between the two gaits cannot occur if the 
system is to remain in these regions. In other words, 

xeEgAyeEZ =► 

(2) 

K a {y)iEl A Wp{x)iEZ Va,/3. 

In this study, we will show how transitions between gaits are found at points outside 
these stable regions. The transitions require the selection of the angle of attack; therefore 
CAAP's are not suitable for this task. We will also show evidence indicating that it is 
possible to find an angle of attack 6 that maps a point into a stable region, e.g. x 
E^ A (3 9, y | y ^ x, y E E^, y = lZg(x)) . Additionally, we introduce the concepts of partial 
stability and viability that will be useful in the construction of the transitions presented 
herein. 

This paper is organized as follows. In section [Til we describe the models used for our 
simulations, their representation in state variables and the definition of the discrete map. 
Next, in section we introduce the new concepts, and we show the regions where the 
transitions between gaits exist. Later, in section HVl we discuss about the requirements of a 
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FIG. 1. (Color online) Illustration of the evolution of the SLIP model for running and walking. 
The mass is represented with a filled circle. The color of the fill indicates touchdown event (black), 
takeoff event (white), and the crossing of the section (pink (grey)). The landing leg is pictured 
with a thick solid line, and the leg at takeoff is represented with a blurred line. Due to the passive 
properties of these models, control is necessary only during the swing of the leg, i.e. during free 
fall while running and from point A to B while walking. 
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controller for the system and the implications for robot design and bipedal locomotion. We 
conclude the paper in section [V] with our conclusion. 



II. METHODS 

As explained previously, we use the SLIP model to study bipedal gaits. We adopt the 
framework in [12J , which is described in the language of hybrid dynamical systems. There- 
fore, we reintroduce some notation and definitions. 

To represent the different phases of a gait, the model is segmented into three sub- models. 
We will call these sub-models charts[4] or phases see FigJU Each chart represents the motion 
of a point mass under the influence of: only gravity (ff-chart or flight phase), gravity and 
a linear spring (s-chart or single stance phase), gravity and two linear springs (d-chart or 
double stance phase). The point mass represents the body of the agent and the massless 
linear springs model the forces from the legs (Fig{T]). A trajectory switches from one chart to 
another when some real valued functions evaluated on it cross zero (event function^, 
We define a running gait as a trajectory that switches from the s-chart to the ff-chart and 
back to the s-chart. A walking gait is defined as a trajectory that switches from the s-chart 
to the d-chart and back again to the s-chart. Switches from the ff-chart to d-chart or vice 
versa are not included in this study. 

A. Equations of motion in each chart 

The motion in all the charts is governed by a system of ordinary differential equations: 



k = F t (a) , (3) 

— * — * 

where X is the vector of state variables and is a force function characteristic of each 
chart. Since all forces are conservative, the energy of the system is constant. For the ff-chart 
the state is described by the Cartesian coordinates of the position of the point mass and its 
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velocity X ff = (x,y,v x ,v y ) T , 



X 



ff 



fv x \ 


\~9j 



(4) 



where g is the acceleration due to gravity. 

The state in the s-chart is represented in polar coordinates X s = (r, 9, r, 9j , where r is 
the length of the spring and 9 is the angle spanned by the leg and the horizontal, growing 
in clockwise direction. Thus, the equations of motion are: 



X s 



r 


^{r -r) + r9 2 - gsin0 
1 ^2f9 + g cos O^j 



(5) 



It is important to note that 9{txD) — a , i- e - the angular state at the time of touchdown is 
equal to the angle of attack. The parameter r defines the natural length of the spring. 

In the d-chart the state is also represented in polar coordinates X d = (r, 9, r, 9^j , with 
the origin of coordinates in the new touchdown point. The motion is described by: 
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(6) 



r & = \Jr 2 + x 2 — 2rx d , cos 9, 



(7) 



where x a is the horizontal distance between the two contact points and is the length of 
the back leg. 



B. Event functions 



Event functions are functions on the phase space of the system. An event occurs when 
the trajectory of the system intersects a level curve of the event function. At the time of the 
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event, the current state of the system is mapped to the state of another chart. Some event 
functions are parameterized with the angle of attack and the natural length of the springs. 
Switches from the ff-chart to the s-chart are defined by: 

/ \ I y — r cos a = 

Fff^s \X ff ,a,r j : < , (8) 

K 

which means that the mass is falling and the leg can be placed at its natural length with 
angle of attack a. Therefore, the motion is now defined in the s-chart. The switch in the 
other directions is simply: 

J r a -+ ff (x„r ):r-r = 0. (9) 



These are the only two event functions involved in the running gait. The map from one 
chart to the other is defined by: 

x = —rcos6 y = rsin9. (10) 

It is important to have in mind that the origin of the s-chart is always at the touchdown 
point. 

For the walking gait, we have to consider switches between single and double stance 
phases. From the s-chart to the d-chart, we have: 

/ _ \ I r sin 9 — r cos a = 
F s ->d [X s ,a,r J : <j , (11) 

l*>i 

which is similar to (jHJ) with the additional condition that the mass is tilted forward. Addi- 
tionally, if we consider the sign of the radial speed, we differentiate between walking gait W 
with f < and Grounded Running gait QTZ, with r > 0. 

The switch from the double stance phase to the single stance phase is defined by: 

J r d +.(Xd,r o y.r l ,-ro = 0, (12) 

with r# as defined in (EJ). The map from the d-chart to the s-chart is the identity. In the 
other direction we have: 

rd = r 9 d = a, (13) 

= r cos a — r s cos 9 S , (14) 



where the subscripts indicate the corresponding chart. 

If the system falls to the ground (y < 0), attempts a forbidden transition (e.g. d-chart to 
ff-chart), or renders v x < (motion to the left, "backwards"), we consider that the system 
fails. 

C. Simulation of the dynamics 

The state of the model is observed when the trajectory of the system intersects the section 
defined by S : 9 = 7r /2. In this way, the map lZ a : S — >■ S transforms points through the 
evolution of the system from the s-chart to the ff-chart and back again to the s-chart using an 
angle of attack a. Similarly, the map W a '■ <S — >■ S transforms points through the evolution 
of the system from the s-chart to the d-chart and back again to the s-chart using an angle 
of attack a. 

All initial conditions are given in the S section and in the s-chart, i.e. only one leg 
touching the ground and oriented vertically. Moreover, all the initial conditions are given at 
the same total energy. The results are visualized using the values of the length of the spring 
r and the radial component of the velocity which, in S, equals the vertical speed f = v y (v x 
is obtained from these values and the equation of constant energy). It is important to note 
that all possible values of r, v y and v x , for a given value of the total energy E, lay on an 
ellipsoid. Besides, there is a transformation that maps the ellipsoid to a sphere. This can 
be shown as follows: the total energy in the section is, 



E = -k (r - r) 2 + -m (vl + Vy) + mgr 



(15) 



Defining the parameters 




(17) 



(16) 



the new variables 
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(18) 
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(19) 




(20) 
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TABLE I. Values used for the simulations presented in this paper. 



Description 


Name 


Value 


Mass 


m 


80 kg 


Elastic constant of linear springs 


k 


15 kNm 


Rest length of linear springs 




1 m 


Total energy 


E 


820 J 


Acceleration due to gravity 


a 


9.81 m/s 2 


Angle of Attack 


a 


from 55°to 90° 



transform equation (TTojl into 



L 2 = vl + v\ + f 2 (21) 



which defines a sphere. Therefore, all initial conditions of f and v y with constant energy, 
are defined inside a circle. A Delaunay triangular mesh was created in the circle with 65896 
initial conditions as vertices (131245 triangles). Each vertex was transformed using TZ a , 
QTZ a and W a with 400 values of a G [55°, 90°]. To compute the evolution of an arbitrary 
initial condition, we used bilinear interpolation in the triangles of the mesh. 

The model implementation and data analysis were carried out in MATLAB(2009, The 
MathWorks), GNU Octave Q and Matplotlib [15j . Simulations were run for constant energy, 
using the step variable integrator ode45 (relative tolerance: 1 x 10~ 6 and absolute tolerance: 
1 x 10~ 8 ). Table fl] shows the values of the parameters used. 



III. RESULTS 

In this section, we present the results of the analysis on the data collected from the 
models as described in section III CI Aiming to define a controller, we introduce some im- 
portant properties of the dynamics of each gait, namely finite stability for a given CAAP 
and viability. 
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A. Finite stability and Viability 



Finite stability describes the set of initial conditions where the system can do a maxi- 
mum amount of steps (sequential applications of the map) before failing, using CAAP. For 
example, we can define for W 

= {x\ x G S A (3a\ y = W n a (x) , n > 1, y G S)}. (22) 

That is, at a given state x = (r, v y ) in S there is a CAAP (a) such that the system can 
do at most n steps before failing. The region E^ are all the points in the section where 
applying W produces a failure. The existence of E^ implies that a controller of the system 
may not need to take a decision at each step. In addition, the controller may exploit this 
alleviation by planning future angles of attack. Viability describes how easy is to choose the 
future angle of attack. The level of ease is measured in terms of the size of the interval of 
angles that can be chosen to avoid a failure of the system. For the running gait this region 
is defined as: 



V R (Aa) ={x\ x eSA 

(23) 

(3a G I a , \\I a \\ > Aa | y = TZ a (x) , y G <S)}, 

where I a denotes a real interval and || ■ || measures its length. In a real system, it is required 
that a viable angle of attack exists for a definite interval, since real sensors and actuators 
have a finite resolution and are affected by noise. 

? jg. [2] shows the finite stability regions for each gait. The stable region of TZ, as reported 



in 



12| (v y = 0) is not visible. Although E^ may have some area of attraction, due to the 



resolution we used for the angles of attack (described in section III Cj) we do not see it in our 



results. Based on results not presented here, we estimate that the resolution in the angle of 
attack to detect such basin for the current energy is ~ 10~ 4 . In despite of the low resolution 
in the angles, the system can perform an average of 10 steps in TZ, and at least 25 steps 
(maximum calculated) in Q1Z and W. This means that running is more difficult at this 
energy level than the other two gaits. Particularly for QTZ and W, we see that there is a 
plateau with the maximum number of steps. This is the evidence of the self-stable regions 
of these gaits, and the plateau is related to the basing of attraction of that region. 
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FIG. 2. (Color online) Finite stability regions. The figures show initial conditions for 1Z, Q7Z and 
W that can do multiple steps under CAAP before failing. A region in white corresponds to Eq for 
gait i. 

Fig. [3] shows the V l (Aa) regions for each gait i. Comparing with Fig. [2], we see that 
in general long partial stability implies wider options for the angle of attack. Particularly, 
transitions are found near these regions of high viability and long partial stability, as will 
be described in the next section. 

Running Ground Running Walking 
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FIG. 3. (Color online) Viability regions for each gait. The figures show the range of angles of 
attack that can be selected in each initial condition that allows the system give at least one more 
step. Colors indicate the size of the window, spanning from 0°to 10°. 

Fig. H] shows one of the strongest results presented here. For each gait i, there is at least 
one angle of attack that maps the current state of the system into E^, and this angle exists 
for an extense region of S. This implies that if we consider control policies with variable 
angle of attack, almost any point in the section can be rendered stable. For this region the 
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optimal control policy requires two angles: the first one maps the point to E^; the second 
angle, keeps the system in this region. 
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FIG. 4. (Color online) Points that can be mapped to stable regions in one step. The figures 
show the initial conditions that can be mapped to a small neighborhood of the stable region E^, 
\v y \ < 1 x 10~ 3 (vy = 0, dashed horizontal lines). Color indicates the angle chosen. Regions V l {2°) 
are marked with solid lines. 



B. Transition regions 

As it was shown in the previous section, the only way of producing transitions between 
gaits is to put the system in a region with finite stability (due to the empty intersection of 



the regions reported in 12j . see FigHj). In Fig. [5] we show transitions starting at E % n and 
arriving at V j (2°) for i ^ j and (i ->■ j) = {(R -» GR), (GR -> W), (W -> GR), (W -> R)}. 
We show the transitions that will be used in the next example, however transitions between 
two any gaits are possible. It shall be noticed that wherever two regions of different gaits 
intersect, the transition is trivial. 

Finally, Fig. [6] and Fig. [7] show one example of three transitions for a given initial condi- 
tion. The trajectory has a total of 26 steps and the angle sequence is 

a = (81.886 5 , 88.500, 62.400, 72.350, 71.100 3 , 71.000, 

(24) 

74.400, 72.130, 74.000 4 , 78.000 2 , 76.500, 69.000, 81.728 4 ) 

where the exponent indicates how many times the angle was used. The path of the center 
of mass in the Cartesian plane is also shown in the figures. 
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FIG. 5. (Color online) Transitions regions landing in Aa > 2°. All the initial conditions that have 
a future inside the region with Aa > 2° of the objective gait are plotted with black dots. The same 
region of the starting gait is given as a reference and appears shaded. Colors in the objective region 
indicate the angle of attack used to perform the transition. Wherever two regions of different gaits 
intersect, the transition is automatically given. 
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FIG. 6. (Color online) Transition sequence. The plot shows a trajectory with three transitions. 
The Regions V 1 (2°) are shown shaded with self-stable regions in dotted line. The arrows indicate 
the order of the sequence and the step number is given. The angle of attack sequence is given in 
(IMD- 
All together we have shown that the SLIP model can be easily controlled to present 
transitions between gaits. To find transitions we must search for an intersection between 
the future of the starting region and the desired objective region. Depending how these 
regions are defined, it may be the case that multiple steps are required to achieve a successful 
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FIG. 7. Transition time series. The figure shows the motion of the point masqj in the plane is 
shown together with the crossing of the section (filled circles [6]). Transition points are indicated 
with a vertical line. 

a An animation of these transitions can be seen in 



htt p : / / www . ifi . uzh. ch/arvo / ailab / people/hamarti / GaitT . avi 



transition. 



IV. DISCUSSION 

There are two important aspects regarding the viability regions. First, it is important 
to notice that V l (Aa) enclose the E^ region, and the points that can be mapped to stable 
regions in one step (Fig. H|) . Second, as it can be seen in Fig. [3l the bigger the range of 
the angle of attack is, the smaller the viability region is. We can take advantage of these 
properties to stabilize the system more easily. The selection of an appropriate Aa e.g. 2° 
defines a set of V l (Aa) inside the section S, where the controller has at least a range of 2° 
to select an appropriate angle of attack. Moreover, the agent can select conservative angles, 
step by step, to bring itself to the region (Fig. [3]). 

Despite the relief to the controller induced by the viability region, the selection of the 
Aa can generate regions that do not intersect; e.g. in Fig. H]we can see that V 1 (2°) does 
not intersect any other region, which makes the gait transition more difficult to carry out. 
In order to cope with this situation, we look at the future of all the initial conditions in E l n . 
As it is presented in Fig. we found that there are some initial conditions, that under a set 
of angles of attack, are mapped from E\ to E 3 n (e.g. E^ to E^ R ). What is also important 
is that the region where we can find these initial conditions are inside the viability region 
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(Fig. ED. 

In these terms, the controller has two purposes. First, based on the state on the S 
section, it has to select the gait, and the angle of attack to keep the agent stable. Thus, the 
controller needs to have the knowledge of all the V R (Aa), and the desired Aa to identify 
which gait has to be selected; the angle of attack can be selected based on the gait model. 
Second, the controller has to be able to produce gait transition when it is needed. Hence, 
the transition regions should be known by the controller and with a model of the gait, the 
angle of attack required can be selected. We expect that this approach can be used to handle 
uneven terrain, given that these irregularities can be modeled (under certain restrictions) as 
a change in energy. 

All these results are conditioned to the selection of the S section. This means that we 
are analyzing the system in only one point in the whole trajectory. From what we see in 
these results, in some regions the trajectories are very close. It would not be a surprise that 
these trajectories of 1Z, W, and QTZ cross each other in another point along their continuous 
evolution, but given that we are looking just at the S section, this cannot be anticipated. 
Nevertheless, the selection of this section establishes the angle of attack as a natural control 
action to stabilize the system and to generate the transitions. 

V. CONCLUSION 

In the present study we have taken advantage of the perspective of hybrid dynamical 
systems to represent locomotion as a process generated by several charts. Although, this 
view makes evident a bigger set of connections among the charts, in this paper we take 
into account a small subset (s-chart to ff-chart, and s-chart to d-chart) which allow us to 
discover new alternatives to perform gait transitions. The development of the maps W a , 
GTZ-a, T^a i s fundamental to identify important regions in the S section that bring the system 
to stable locomotion and to a gait transition. The present results bring new ideas about 
plausible mechanisms that biped creatures could use to carry out gait transitions and stable 
locomotion. These mechanisms exploit the passive dynamics of the system, which reduces 
the amount of energy needed to control the system. These features are also present in biped 
machines with compliant legs, and as suggested in this paper, these mechanisms can be 
exploited to develop stable gaits and gait transitions. 
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